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Results on heat current, entropy production rate and entanglement are reported for a quantum 
system coupled to two different temperature heat reservoirs. By applying a temperature gradient, 
different quantum states can be found with exactly the same amount of entanglement but different 
purity degrees and heat currents. Furthermore, a nonequilibrium enhancement-suppression transi- 
tion behavior of the entanglement is identified. 
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I. INTRODUCTION 

Quantum thermodynamics 0, Q is starting to throw 
light on universal behaviors of nanosystems. Specifi- 
cally, new possibilities arising in nonequilibrium situa- 
tions, with dominant quantum coherences, are emerging 
Currently, non-local quantum correlations (entan- 
glement) are being considered in a vast variety of sce- 
narios as they are key ingredients for novel and non- 
conventional forms of communication, information pro- 
cessing and computation. For potential large-scale ap- 
plications, where condensed matter systems are of prime 
importance, thermal interactions with specific environ- 
ments are unavoidable. Thus a clear connection between 
quantum information aspects and thermal magnitudes 
has to be elucidated. Furthermore, the biological frontier 
of physics imposes to address the question of quantum 
features survival in noisy as well as in nonequilibrium 
conditions [3]. 

For quantum systems in contact with heat reservoirs at 
unique and fixed temperature, the equilibrium thermal 
entanglement has been extensively studied 0, 0, S 0- 
Until now, most of the emphasis in the study of ther- 
mal entanglement has been confined to equilibrium situa- 
tions. Entanglement in nonequilibrium quantum systems 
has been scarcely considered and thus a proper descrip- 
tion of thermal entanglement in the presence of matter 
and/or energy currents is still lacking. Recently, Eisler et 
al. [9( has calculated the Von Neumann entropy of a block 
of spins in a XX spin chain in the presence of an energy 
current, showing that an enhancement of the amount of 
entanglement due to an energy current is possible. The 
energy current is modelled by adding an extra term to 
the spin chain Hamiltonian for simulating a steady-state 
current in a thermodynamic closed system. However, the 
issue of entanglement behavior in true nonequilibrium 
conditions of a thermodynamic open system remains un- 
touched. 

Two coupled qubits in thermal contact with different 
heat baths is a system not only of theoretical interest but 
a common place in nanophysics. In semiconductor quan- 



tum dots the transfer of quantum information between 
nuclear spins and electronic spins has been recently con- 
sidered 13, 11 , 12 1 . The nuclear spin is generally weakly 
coupled to its environment while the electronic spin is 
strongly coupled to a great variety of degrees of freedom 
within the solid. In this way, the effective environments 
are different for both kind of spins. Besides that, nu- 
clear magnetic resonance techniques allow the cooling of 
nuclear spins in a controlled manner (l3l . [3] without sig- 
nificatively affecting the electronic spins, thus creating 
two reservoirs at effective different temperatures. On the 
other hand, superconductor qubits can be easily designed 
to be coupled to different environments. For instance, in- 
ductively coupled superconductor flux qubits in contact 
with two different environments has been recently ana- 
lyzed in Ref . [15[ . The aim of the present paper is to cor- 
relate thermodynamical nonequilibrium steady-state fea- 
tures with entanglement properties of quantum nanosys- 
tems. In doing so, we shall consider a quantum system 
in a nonequilibrium condition for which the amount of 
entanglement can be exactly evaluated: two interacting 
qubits (spins) in contact with two heat reservoirs at dif- 
ferent temperatures. In this case, entanglement can be 
evaluated for any mixed state by using the concurrence 
[lij ]. Indeed, the model system to be considered in the 
present work should be useful for a large variety of phys- 
ical set ups aiming to explore the relationship between 
quantum informational entropy and thermodynamic en- 
tropy at the atomic scale. Whether one can reveal uni- 
versal features in irreversible processes of open quantum 
systems is of great significance. 



II. FORMALISM 



The central quantum nanosystem is described by a 
Hamiltonian Q, interacting with two heat reservoirs 
which are assumed to be in a permanent thermodynam- 
ical equilibrium at /3j = l/kgTi,! = 1,2 (/eg = % = 1) 
with internal Hamiltonians i? ? . The total Hamiltonian is 
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then 



H = Q + R 1 +R 2 + S 1 + S 2 



(1) 



where the nanosystem is simultaneously coupled with 
both reservoirs through terms Si y S 2 - The nanosystem 
+ reservoirs is described by a density operator satisfy- 



ing the Liouville equation 
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[H, 7]. We assume 



that the coupling strengths of the central quantum sys- 
tem to the reservoirs are weak so that the full density 
operator 7 can be expressed as ^(t) = p{t)pip 2 where 
each reservoir is described by its own canonical equilib- 



rium density operator pi 



/T rRi {e-^ Ri }} and 



p(t) is the reduced density operator for the quantum sys- 
tem of interest. The couplings nanosystem-reservoirs are 
written as 



& + & = E E Wai. = EE V 1JL 



(2) 



i=i m 



where the nanosystem operators Vj. M are taken to sat- 
isfy [Q, Vj.fj] — ^j, fj.Vj.fi and the operators fj jfl act on 
the reservoir degrees of freedom (j = 1,2). Within the 
framework of the Born-Markov approximation 17|, the 
equation of motion for p(t) is given by 



dp 



i[4p(*)]-EE J ^K-) 

J'=l Vl" 

{[Vj#, [Vl, p]] - (1 - e^)[V^, Vlj}} 



' 2,111 rjj r " / L 'j, 

where the spectral density of the j'-th reservoir is 



(3) 



with /L( r ) = ' 



"/,•,,. {/-;./•;;„-)./•,,,} (4) 
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We will here be concerned with the simplest possible 
scenario where clear relations between informational and 
thermodynamic entropies could be found. To set up our 
model system in a general context, we consider a nanosys- 
tem composed of two interacting qubits as described by 
the Hamiltonian 
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Q = E Y a °' 

a=l 



-) + K'a hz a 2 , z (5) 



where a a . z and crjj denote Pauli matrices. The inter- 
qubit coupling is ferromagnetic when A", A' < and an- 
tiferromagnetic when A, A' > 0. This type of Hamil- 
tonian encompasses three well-known spin models: it 
turns into the isotropic Heiscnberg-like coupling for A — 
A', the isotropic XX-like model for K' = and the 
Ising-like model for K = 0. The eigenenergies and 
eigenstates corresponding to Eq.© are: \s\) = |0, 0) 
( El = -a±a + A"), \s 2 ) = |1, 1) (A 2 = a±a + A"), 
I s 3 ) = cos(0/2)|l,Q) + sin(0/2)|O,l) (E 3 =a-K') and 
|s 4 ) = -sin(0/2)|l,O) + cos(0/2)|O,l) {E A = -a 

with a = \ A 2 



and tan0 = 2A/(ei - e 2 ) 



K 1 ), 
We 



consider each qubit in contact with its own boson heat 
reservoir through a term of the form 



^ = °t E ffn )fi n,j + ^ E 5 « 



(j)*^t 



J = 1,2 (6) 



where ■ creates an excitation in mode n of reservoir j 
with a coupling strength g„ ■ 

The nonequilibrium steady-state density matrix (de- 
signed simply as p from now on), must satisfy -£ = 
-i[Q,p] = in Eq.©, which yields to d{p) + £ 2 (p) = 
where the Lindblad or relaxation super-operators are 
given by 



fL=l 

4 

- E jU) (-^{-VLpV^ + pV^Vl + e-ft"" (Vj^/S ~ ^>^)> ( 7 ) 



for j = 1,2. In the latter expression u>x = E 2 — 
E 3 ,V s ,i = (<5 i)2 cos(0/2) + 5 jA sm{6/2))\s 2 )(s 3 \; u 2 = 
E 2 - Ei,V j<2 = (-<5 i)2 sin(0/2) + <5 ia cos(^/2))| s 2 ) (s 4 1 ; 
w 3 = /•• ) /•■ 1 • I = (* A icos(0/2) + a A2 sin(0/2))|« 3 }<«i|; 
w 4 = - (-5 jt ism(6/2)+5 j , 2 cos(6/2))\s 4 ,)(s 1 \ 

and jU^—Ufj) = e^ w " J^(w M ). Two limiting cases 
can be easily analyzed: (i) No inter-qubit coupling, 



K = K' = (0 = 0) which leads to = w 4 = e 2 and 
cj 2 = u>3 = e\. Each qubit reaches a local equilibrium 
with its own heat reservoir yielding to a direct product 
form of the density matrix and thus no- entanglement, (ii) 
Coupled qubits, A, K' 7^ 0, in contact with two indepen- 
dent reservoirs at identical temperatures, (3\ = f3 2 = f3. 
A reduced density matrix results which has the thermo- 
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dynamical canonical form for a system described by in- 
ternal Hamiltonian Q at equilibrium with a thermal bath 
at inverse temperature /3, as it should be. 

III. RESULTS AND DISCUSSION 

Consistently with the Born-Markov approximation, 
we adopt a Weisskopf-Wigner-like expression such as 
J^\lo^) = Tj(ujfj,)rij(ujfj,) where Tj^^) depends on both 
the nanosystem-jth-reservoir coupling strength and the 
reservoir internal structure. On the other hand, nj(u>fi) = 
{eP^f — denotes the thermal mean value of the num- 
ber of excitations in reservoir j at frequency uj^. For the 
sake of simplicity, we take identical and frequency inde- 
pendent couplings, thus Tx(u>) = T 2 (u}) = T. 

From Eq.© the nonequilibrium steady-state density 
matrix is obtained as given by the diagonal matrix p = 
rfia<7{pi ; i, p2,2, P3,3, hi the basis of eigenstates of Q. 
Although, it can be analytically expressed we will not go 
here into the details as its explicit form is cumbersome 
[lij ]. Instead, we shall analyze some important special 
situations. 



A. Symmetric case, ei — £2 = e 

In this case p can be written in terms of a simple uni- 
versal function e(uo) = "iM+"2M <- ^ (energies and 
temperatures in units of interqubit coupling K = 1). In 
the strong coupling case (e < 1) we found 

ei / e 2 \ / ei\ e 2 

where ej = e(ujj) with uj\ = W4 = |e — 1| and w 2 = 
u>3 = e + 1. In the weak coupling case (e > 1) the fol- 
lowing interchanges have to be made: pi.i <-> ^4,4 and 
P2.2 P3.3- Thus, the nonequilibrium concurrence is 
C = 2Max{0, 1/93,3 - P4,4|/2 - y/pTJp^}. 

Let us first discuss the equilibrium (T± = T2 = T) 
thermal entanglement behavior for a system governed 
by Hamiltonian ((5]) 0, 0, S S]- An analytical expres- 
sion can be found for the equilibrium concurrence as 

C eg(T) = 2cosh( "°/ 2 ( rSh(l a /2T) ■ This last expression 
is interesting because it implies an universal form (in- 
dependent of e) for the sudden death of the equilibrium 
concurrence at the temperature T c = 1.1346. For e < 1, 
the two-qubit concurrence decreases from 1 to as the 
temperature increases up to T c ; for e > 1, the concurrence 
increases from to some maximum before vanishing at 
T c . It is also known that the concurrence decreases mono- 
tonically as the qubit splitting increases for any temper- 
ature and vanishes exponentially with increasing e. All 




FIG. 1: (Color online) Quantum state parameter space (ei, 
ea). The black line denotes equilibrium states. Green line: 
entangled-unentangled border for e < 1. Red line: entangled- 
unentangled border for e > 1. Circles represent the state shift 
under a temperature gradient from the same mean tempera- 
ture Tm = L Blue circles: e = 1/3, brown circles: e = 3. 

these features are independent of the specific nature of 
the reservoirs. 

The behavior of quantum states, for any qubit internal 
splitting and reservoir temperatures, can be displayed in 
a single plot e±-e2, as it is shown in Fig. 1. Although, 
the quantum state variation is described by the same 
curve for e and 1/e, the border lines separating the en- 
tangled from the unentangled regions are different: the 
green curve corresponds to e < 1 while the red curve cor- 
responds to e > 1. The black line represents equilibrium 
quantum states for both e = 1/3 as well as for e = 3. 
The shifting of the quantum state with the temperature 
gradient, AT = T\ — T 2 , is depicted by the circles, for 
the same average temperature, Tm — iT\ +T 2 )/2 = 1. It 
is evident that the temperature gradient shifts the state 
from the entangled zone to the unentangled zone. How- 
ever, this behavior can be reversed at low temperatures 
for e > 1 as it is to be discussed below. 

In the linear nonequilibrium limit (LNEL), AT <C 
1, the concurrence can be written as C(AT) = 
Max{0, C eg (Tj^f) — qAT 2 } where the coefficient a is 
a function of the average temperature as well as the 
qubit internal splitting. The equilibrium concurrence 
is displayed in the insets of Figs. 2-a and 2-b. In the 
strong coupling limit (e < 1) a > thus the concur- 
rence is always a decreasing function of the temperature 
gradient AT. By contrast, in the weak coupling limit 
(e > 1) there is a transition mean temperature for which 
a changes the sign. Thus, a low temperature region can 
be found where a < for which a gradient tempera- 
ture produces an increasing of the concurrence as com- 
pared with the equilibrium case. The degree of mixing 
of the quantum state can be characterized by the linear 
entropy as defined by S L = (4/3) (1 - Tr{p 2 }). The low 
AT limit of the linear entropy can also be expanded as 
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FIG. 2: (Color online) LNEL coefficients a and a', denot- 
ing the second order in AT variation of the concurrence and 
linear entropy, respectively: (a) e = 3, (b) e = 1/3. Insets: 
Equilibrium concurrence. 

S L (AT) = S L ,e g (T M ) + a' AT 2 + 0(AT 4 ). The coeffi- 
cient a' > 0, for both interqubit coupling cases, is also 
illustrated in Figs. 2-a and 2-b. Note that while a tem- 
perature gradient can produce, in a limited temperature 
interval, an enhancement of the concurrence it always 
yields to a more mixed state. This result will permit to 
prepare a great variety of quantum states with practi- 
cally any combination of entanglement and purity degree 
by varying the temperature of only one heat reservoir. 

The relationship between a nonequilibrium thermody- 
namical quantity such as the heat current and a central 
quantum information concept such as entanglement is 
now addressed. We start by calculating the heat current 
as Jj = Tr{QCj(p)} which in the symmetric case 
yields to 

Ji = J[wi(l-ei)(na(wi)-ni(a;i))-(1^2)] (9) 

and J 2 = -Ji = -J. In LNEL, AT < 1, the 
Fourier's law is well verified, i.e. J — kAT, with k 
the thermal conductance depending on the qubit inter- 
nal splitting and mean temperature [lij ]. The evolution 
of heat current and concurrence, as AT increases is il- 
lustrated in Fig. 3 (obviously J = for AT = 0). 
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FIG. 3: (Color online) Heat current J and concurrence evo- 
lution for different mean temperatures Tm- Inset: Linear 
entropy and concurrence. Each point corresponds to a AT 
value. 
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FIG. 4: (Color online) Thermodynamical entropy production 
rate and concurrence evolutions for different mean tempera- 
tures Tm. Each point corresponds to a AT value. 

Clearly, for temperatures for which a < 0, see Fig. 2- 
a, an enhancement of the concurrence is possible by 
applying a temperature gradient. Based on this gen- 
eral concurrence's behavior, we conclude that in LNEL, 
J {AT) = -^\C eq - C(AT)| 1 /2 as it i s clearly observed 
in Fig. 3. A remarkable point to be noted is the possibil- 
ity of constructing nonequilibrium quantum states with 
identical concurrence, as that for the equilibrium case, 
but carrying a heat current. It is also evident from Fig. 
3 that the relation between heat current and concurrence 
becomes independent of Tm as the temperature gradient 
increases. The correlation between quantum linear en- 
tropy S l and concurrence is also shown in the inset of 
Fig. 3, confirming the fact that a gradient temperature 
will always increase the mixing degree of the quantum 
state. Although, the amount of entanglement is small in 
those cases, it can be significantly increased by distilla- 
tion protocols. 

Any heat current produces an amount of thermody- 
namical entropy proportional to the heat which is car- 
ried on and inversely proportional to the temperature 
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FIG. 5: (Color online) Heat current J and concurrence evolu- 
tion for different mean temperatures Tm in the non-symmetric 
case, ei = 8 and £2 = 3. Inset: ei = 3 and 62 = 8. Each point 
corresponds to a AT value. 

of the reservoir from which the heat is extracted (or 
injected). Thus, the thermodynamic entropy produc- 
tion rate in our system can be written as dS/dt = 
Ji(Tf 1 - T^ 1 ) [19]. In LNEL, the rate of entropy pro- 
duction shows a linear dependence with the concurrence 
as §■ = -M-\C e „ - C(AT)\ as it is shown in Fig. 4. 

Q M 

The entropy production rate, like the heat current and 
linear entropy, is always different for two different val- 
ues of AT corresponding however to the same amount of 
entanglement. 

B. Non-symmetric case, ei > C2 

We first consider the high temperature reservoir (Ti) 
is in direct contact with the large splitting qubit, ei, 
and the low temperature reservoir (T 2 ) is in contact 
with the small splitting qubit, e 2 . Modifications to 
equilibrium values of physical magnitudes such as the 
concurrence and the linear entropy are now of first or- 
der in AT instead of order AT 2 , as it was the case 
for the symmetric set up. This implies that in LNEL 
J (AT) ~ \C eq - C(AT)\, as is illustrated in Fig. 5. 
Note that at low temperature it arises the possibility of 
finding up to three quantum states with the same con- 
currence but carrying different heat currents. Switching 
to the inverse connection between qubits and reservoirs 
(Ti «-> T 2 ), the heat current dependence on the concur- 
rence is completely modified. In this latter case, high 
temperature bath in contact with the low splitting qubit, 
the heat current is substantially decreased but the con- 
currence can be enhanced by the temperature gradient, 
as it is shown in the inset of Fig. 5. We conclude that a 
qubit splitting asymmetry brings an interesting new con- 



trol parameter for engineering nonequilibrium thermal 
quantum states. 

IV. CONCLUSIONS 

In summary, we have demonstrated that under 
nonequilibrium thermal conditions a versatile scenario 
for tailoring heat carrying quantum states with a well 
specified amount of entanglement is feasible. A temper- 
ature gradient has been shown to produce increasing or 
decreasing entanglement depending on the internal cou- 
pling strength within a nanosystem. Physical realizations 
of the model system we addressed are provided by a large 
number of physical systems such as nuclear spins in quan- 
tum dots and superconducting qubits. Therefore, the re- 
sulting insights can serve as useful recipes for realistic 
quantum information processors in noisy and nonequi- 
librium environments. 
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